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ABSTRACT 



^ ' Context. Precise measurements of surface abundances of extremely low metallicity stars have recently been obtained and provide 
new constraints for the stellar evolution models. 

Aims. Compute stellar evolution models in order to explain the surface abundances observed, in particular of nitrogen. 
Methods. Two series of models were computed. The first series consists of 20 Mq models with varying initial metallicity 
(Z = 0.02 down to Z = 10~ 8 ) and rotation (v mi = 0- 600 kms" 1 ). The second one consists of models with an initial metallicity 
of Z = 10~ 8 , masses between 9 and 85 Mq and fast initial rotation velocities (tij n j = 600 — 800 kms -1 ). 

Results. The most interesting models are the models with Z — 10~ 8 ([Fe/H]~ —6.6). In the course of helium burning, carbon 
and oxygen are mixed into the hydrogen burning shell. This boosts the importance of the shell and causes a reduction of the 
CO core mass. Later in the evolution, the hydrogen shell deepens and produces large amount of primary nitrogen. For the most 
massive models (M > 60 Mq), significant mass loss occurs during the red supergiant stage. This mass loss is due to the surface 
Q ■ enrichment in CNO elements via rotational and convective mixing. The 85 Mq model ends up as a WO type wolf-Rayet star. 

' Therefore the models predict SNe of type Ic and possibly long and soft GRBs at very low metallicities. 
C/} The rotating 20 Mq models can best reproduce the observed CNO abundances at the surface of extremely metal poor (EMP) 
C3 • stars and the metallicity trends when their angular momentum content is the same as at solar metallicity (and therefore have 
| an increasing surface velocity with decreasing metallicity). The wind of the massive star models can also reproduce the CNO 
abundances of the most metal-poor carbon-rich star known to date, HE1327-2326. 
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1. Introduction 

Precise measurements of surface abundances of extremely 
metal poor (EMP) stars have recentl y be en obtained by 
ICavrel et afll^OO^ . Ispite et alJ (|2005^ and llsraelian et all 
ll2004ah . These pro vide new constraints for the stellar evo- 
lution models (see IChiappini et all l2005t iFrancois et alJ 
l2004t IPrant^oll200a) 7 The most striking constraint is the 
need for primary 14 N production in very low metallic- 
ity massive stars. Other constraints arc an upturn of the 
[C/O] ratio with a [C/Fe] ratio (see for example Fig. 14 
for C/O and Fig. 13 for C/Fe in lSpite et ailfcOOSl) about 
constant or slightly decreasing (with increasing metal- 
licity) at very low metallicities, which requires an in- 
crease (with increasing metallicity) of oxygen yields be- 
low [Fe/H]^ -3. About one quarter of EMP stars are car- 
bon rich (C-rich EMP, CEMP stars). iRvan et all l|2005(l 
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and lBeers k, Christliebl l|2005|) propose a classification for 
these CEMP stars. They find two categories: about three 
quarter are main s-process enriched (Ba-rich) CEMP 
stars and one quarter are enriched with a weak compo- 
nent of s-process (Ba- normal). The t wo most metal poor 
stars known to date, HE1327-2326 llFrebel et all 120051 
Aoki i et all l200fil and HE 01 07-524o" lchristlieh et ail 



2004) are both CEMP stars. These stars are believed to 



have been enriched by only one to several stars and the 
yields of the models can therefore be compared to their ob- 
served abundances without the filter of a galactic chemical 
evolution model (GCE). 

The evolution of low metallicity or m etal free stars 
is not a new subject (see for example IChiosil Il983t 
lEl Eid et al.ll983HCarr et al.ll984tlArnettll996() ."The ob- 
servations cited above have however greatly increased the 
interest in very metal poor stars. There are many recent 
works st udying the evolution of metal free (or almost) 
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massive ( Heeer fc Wooslevll2002t iLimongi fc ChieffilEoOE 
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lUmeda fc Nomotol 120051 iMevnet et al.ll2006l). intermedi- 
ate mass (ISiess et al ll2002tlHerwidl2004tlSuda et al 1 1 2004 ; 



Gil-Pons et al.ll2005tl and low mass l|Picardi et all 12004 ; 
Weiss et al.ll2004 stars in an attempt to explain the ori- 
gin of the surface abundances observed. In this work pre- 
supernova evolution models of rotating single stars were 
computed with metallicities ranging from solar metallicity 
down to Z = 10~ 8 to study the impact of rotation and to 
see which initial rotation velocities can lead to the chem- 
ical enrichments observed. In Sect. 2, the model physical 
ingredients and the calculations are presented. In Sect. 
3, the evolution of the models is described. In Sect. 4, the 
stellar yields of light elements are presented and compared 
to observations. In Sect. 5, the conclusions are given. 

2. Description of the stellar models 

The stellar evolution model used to calculate the ste l- 
lar models is described in detail in Irlirschi et alJ l|2004|) . 
Convective stability is determined by the Schwarzschild 
criterion. Convection is treated no longer as an instanta- 
neous mixing but as a diffusive process from oxygen burn- 
ing onwards. The overshooting parameter is 0.1 Hp for H 
and He-burning cores and o therwise. The reacti on rates 
are taken from the NACRE I Angulo et al.lll999^ compi- 
lation for the experimental rates and from the NACRE 
website (http://pntpm.ulb.ac.be/nacre.htm) for the theo- 
retical ones. 



2.1. Initial composition 

The initial composition of the models is given in Table 
[U The solar metal licity composition is described in 
Hirs chi et alJ (|2005bfl and only low metallicity composi- 
tions are presented here. For a given metallicity Z (in 
mass fraction) , the initial helium mass fraction Y is given 
by the relation Y = Y p + AY/AZ ■ Z, where Y p is the pri- 
mordial helium abundance and AY/AZ the slope of the 
helium-to-metal enrichment law. Y p = 0.24 and AY/AZ 
= 2.5 were used accor ding to recent determinations (see 
llzotov fc Thuanll2004i for example). For the mixture of 
the heavy elements, the same mixture as the one used to 
compute the opacity tables for Weiss 95 's alpha-enriched 
composition ( Iglesias fc R,ogers1ll99fi|) was adopted. 



Table 1. Initial abundance (in mass fraction) of the chem- 
ical elements followed in the calculations. 



2.2. Mass loss 

Since mass loss rates are a key ingredient for the evolution 
of massive stars, the prescriptions used are summarised 
here. The changes of the mass loss rate s, M, with rotation 
are tak en into account as explained in lMaeder & Mevnetl 
(2000a). As reference mas s loss rates, the adopted mass 
loss rates are the ones of IVink et al.1 l)2000t l200l|) . who 
account for the occurrence of bi-stability limits, which 
change the wind properties and mass loss rates. For the 
domain n ot covered by these a uthors the empirical law 
devised bv lde Jager et al.l l)l988j) was used. Note that this 
empirical law, which presents a discontinuity in the mass 



Element 


Z = 1C 




Z = 10" b 


Z = 10 _s 


1 T T 


7.5650D- 


01 


0.759965 


0.759999965 


3 He 


2.5702D- 


05 


2.5440D-05 


2.5437D-05 


4 He 


2.4247D- 


01 


0.23999956 


0.239974588 


12 C 


7.5542D- 


05 


7.5542D-07 


7.5542D-10 


13 C 


9.0930D- 


07 


9.0930D-09 


9.0930D-12 


14 N 


2.3358D- 


05 


2.3358D-07 


2.3358D-10 


15 N 


9.2242D- 


08 


9.2242D-10 


9.2242D-13 


16q 


6.7105D- 


04 


6.7105D-06 


6.7105D-09 


17 


2.7196D- 


07 


2.7196D-09 


2.7196D-12 


18q 


1.5162D- 


06 


1.5162D-08 


1.5162D-11 


20 Ne 


7.8366D- 


05 


7.8366D-07 


7.8366D-10 


22 Ne 


6.3035D- 


06 


6.3035D-08 


6.3035D-11 


24 Mg 


3.2475D- 


05 


3.2475D-07 


3.2475D-10 


25 Mg 


4.2685D- 


06 


4.2685D-08 


4.2685D-11 


26 Mg 


4.8956D- 


06 


4.8956D-08 


4.8956D-11 


28 Si 


3.2769D- 


05 


3.2769D-07 


3.2769D-10 


32 g 


1.8897D- 


05 


1.8897D-07 


1.8897D-10 


36 Ar 


1.9797D- 


06 


1.9797D-08 


1.9797D-11 


40 Ca 


5.1728D- 


06 


5.1728D-08 


5.1728D-11 


44 Ti 













48 Cr 













52 Fe 













56 Ni 














flux near the Humphreys-Davidson limit, implicitly ac- 
counts for the mass loss rates of LBV stars. For the non- 
rotating models, since the empirical values for the mass 
loss rates are based on stars covering the whole range 
of rotational velocities, one must apply a reduction fac- 
tor to the empirical rates to make them correspond to 
the non-ro tating case. The same r eduction factor was 
used as in iMaeder &; Mevnetl feOOll). During the Wolf- 
Rayet phase the mass loss rates bv lNugis fc Lamersl fcOOCft 
were used. The mass loss rates depend on metallicity as 
M ~ (Z/Zq) - 5 , where Z is the mass fraction of heavy 
elements at the surface of the star. 

The mass loss rates (and opacity) are rather well de- 
termined for chemical compositions which are similar to 
solar composition or similar to a fraction of the solar com- 
position (or of the alpha-enriched mixing). However, very 
little was known about the mass loss of very low metallic- 
ity stars with a strong enri c hmen t in CNO elements until 
recently. IVink fc de Koterl l)2005[) study the case of WR 
stars and find a clear dependence with iron group mass 
fract ions. For red su pergiant stars (RSG), recent studies 
(see Ivan Loonl 120051 and references therein) show that 
dust-driven winds at cool temperature show no metallicity 
dependence for Z > 0.1 Zq. As we shall see later, due to 
rotational and convective mixing, the surface of the star is 
strongly enriched in CNO elements during t he RSG stage- 
It is implicitly assumed in this work (as in lMevnet et alJ 
2006) that CNO elements have a significant contribution 
to opacity and mass loss rates. The mass loss rates used 
depend on metallicity as M ~ (Z/Z & ) - 5 , where Z is the 
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mass fraction of heavy elements at the surface of the star, 
also when the iron group elements content is much lower 
than the CNO elements content. The largest mass losses in 
the present calculations occur during the RSG stage. This 
means that if the independence o f the mass loss r ates on 
the metallicity in the RSG stage ijvan Loonll2005() is con- 
firmed at very low metallicities, the mass loss rate used in 
this work possibly underestimate the real mass loss rate. 
This point surely deserves to be studied in more detail in 
the future. 

A sp ecific treatment for m ass loss was applied at break- 
up fsee iMevnet et al.ll2006|) . At break-up, the mass loss 
rate adjusts itself in such a way that an equilibrium is 
reached between the two following opposite effects. 1) 
The radial inflation due to evolution, combined with the 
growth of the surface velocity due to the internal coupling 
by meridional circulation, brings the star to break-up, and 
thus some amount of mass at the surface is no longer 
bound to the star. 2) By removing the most external lay- 
ers, mass loss brings the stellar surface down to a level in 
the star that is no longer critical. Thus, at break-up, one 
should adapt the mass loss rates, in order to maintain the 
surface layers at the break-up limit. In practise, however, 
since the critical limit contains mathematical singulari- 
ties, it was considered that during the break-up phase, the 
mass loss rates should be such that the model stays near 
a constant fraction (around 0. 95) of the limit. Note that 
wind anisotropy (described in iMaeder fe Mevnetl EoOOa) 
was not taken into account for the present work. 

2.3. Chemical elements and angular momentum 
transport 

The instabilities induced by rotation taken into account 
in this work are meridional circulation and secular and 
dynamical shears. Meridional circulation is an advective 
process and shears are diffusive processes. The equations 
for the transport of angular momentum and chemical ele - 
ments are given in Sect. 2.3 o f lMaeder fc Mevnetl l^nOObft . 
For more details on the equation for the transport of an- 
gular momentum, the reader can refer to lMaeder fc Zahnl 
(1998). The equation for the change of the mass fraction 
of chemical species i is the following: 



change in chemical abundances can be approximated by a 
diffusion process. The coefficient D e g (accounting for the 
combined effect of advection and horizontal turbulence) is 
calculated in this work using the following formula: 



rU(r) 
30D h 



(2) 



dXj 
dt 



Mr 



d 



dM r 



{Attt 2 P ) 2 D 



dM r 



dt 



The second term on the right accounts for composition 
changes due to nuclear reactions. Note that the nuclear 
burning term and the diffusive term are treated in a serial 
way in the simulations. The coefficient D is the sum of the 
different diffusion coefficients: D = D conv + D c ff + D s[ical , 
where D conv is the convective diffusion coefficient, D c s ac- 
counts for the combined effect of advection and horizon- 
tal turbulence a nd J> R h B a r rep resents the sum of secular 
(s ee Eqs. 5.32 inlMaederlll997tl and dynamical (see Eq. 5 
in lHirschi et al.ll2004j) shear coefficients. Although merid- 
ional circulation is an advective process, its effect on the 



where Dh is the coefficient of horizontal turbulence, for 
which the estimate is D h = \rU(r)\ l|Zahnlll992h . This 
equation expresses that the vertical advection of chemi- 
cal elements is severely inhibited by the strong horizontal 
turbulence characterised by Dh- 

Even though there is no free parameter in the prescrip- 
tions above to increase or decrease the importance of the 
coefficients, different authors use differ ent prescriptions fo r 
the various processes (see for example iHesrer et alJl200(i|) . 
The coefficient of horizontal turbulence was also recently 
revised IIMaederl 120031) . The new coefficient was used in 
iMevnet et alJ l)2006j) but not in this work. The impact of 
the new coefficient is discussed briefly in Sect. 13.41 

2.4. Initial rotation 

The value of 300 kms -1 as the initial rotation velocity 
at solar metallicity corresponds to an average velocity of 
about 220 kms -1 on the Main Sequence (MS) which is 
very close to the average observed value (see for instance 
lFukudal F982V It is unfortunately not possible to measure 
the rotational velocity of very low metallicity massive stars 
since they all died a long time ago. Higher observed ratio 
of Be to B stars in the Mage llanic clouds compared to our 
Galaxy I Maeder et al. 1 119991 could point out to the fact 
the stars rotate faster at lower metallicities. Also a low 
metallicity star containing the same angular momentum 
as a solar metallicity star has a higher surface rotation 
velocity due to its smaller radius (one quarter of Zq ra- 
dius for 20 Mq stars). Since there is however not yet firm 
evidence for fast surface rotation velocities at low metallic- 
ities, we explore in this work with 20 Mq models different 
velocities ranging between no rotation and surface veloc- 
ities corresponding to the same total angular momentum 
as in solar metallicity stars. 

In order to compare the models at different metallici- 
ties and with different initial masses with another quantity 
-j^than the surface velocity, the ratio Vi n i/v CL n is used (see 
Table |2J. The critical velocity is reached when the grav- 
itational acceleration is balanced by radiative and cen- 
trifugal forces. The critical velocity is given by the fol- 
lowing formula if the star is far from its Eddington limit 
(r max < 0.64): 



^crit,l = f2 it 



eb 



2GM\ 5 



3 i? p b 



(3) 



i? e b and i?pb are respectively the equatorial and polar ra- 
dius at the break-up velocity. If the star gets closer to the 
Eddington limit, then the following critical velocity has to 
be used: 
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i -r„ 



J crit,2 



4 U ' 



crit,l 



V'(cj) 



MM 

i? 2 



(4) 



where R c (w) is the equatorial radius for a given value of 
the rotation parameter u>. More details can be found in 
iMaeder fc Mevnetl l|2000a[) . Towards lower metallicities, 
Wini/^crit increases only as r -1 / 2 for models with the same 
angular momentum (J), whereas the surface rotational ve- 
locity increases as r _1 ( J ~ vr). The angular momentum 
varies significantly for models with different initial masses. 
Finally, Umi/vcrit is a good indicator for the impact of ro- 
tation on mass loss. 

In the first series of models, the aim is to scan the 
parameter space of rotation and metallicity with 20 Mq 
models since a 20 M p> star is not fa r from the average 
massive star when a ISalneterl l)l955j) like IMF is used. 
For this series, on top of non-rotating models, two initial 
rotational velocities were used at very low metallicities. 
The first velocity is the same as at solar metallicity, i. 
e. 300 km s -1 . The ratio Vi n i/v CI i t decreases with metal- 
licity (see Table El for the initial velocity of 300kms _1 . 
The second u ini is SOOkms" 1 at Z=10" 5 ([Fe/H]~- 3.6) 
and GOOkms" 1 at Z=10~ 8 ([Fe/H]~- 6.6). These values 
have ratios of the initial velocity to the break-up velocity, 
^ini/wcrit around 0.55, which is only slightly larger than 
the solar metallicity value (0.44). The 20 Mq model at 
Z=10 -8 and with 600 km s" 1 has a total initial angular 
momentum Jtot = 3.3 1 52 ergs, which is the same as 
for the solar metallicity 20 Mq model with 300 km s -1 
(J tot = 3.6 10 52 ergs). So even though a star at Z=10~ 8 
with a velocity of 600kms _1 at first glance seems to be 
an extremely fast rotator, it is in fact similar to a so- 
lar metallicity star in terms of angular momentum and 
ratio Uini/f rrit.- In the second se ries of models, following 
the work of IMevnet et alJ l|2006|) . models were computed 
at Z=10 -8 with initial masses of 40, 60 and 85 Mq and 
initial rotational velocities of 700, 800 and 800 km s -1 re- 
spectively. Note that, for these models as well, the initial 
total angular momentum is similar to the one contained 
in solar metallicity models with rotational velocities equal 
to 300 km s -1 . 

3. Evolution 

The evolution of the models was in general followed until 
core Si-burning. The 60 AIq model was evolved until neon 
burning and the 9 Mq model until carbon burning. The 
main characteristics of the models are presented in Table 
□ 



3.1. Metallicity effects 

The effects of metallicity on stellar evolution have al- 
ready been discu s sed in the literature (see for example 
Heger et al.ll2003l IChieffi fc Limongill2004t IMevnet et al. 



1994). But before looking at the very low metallicity mod- 



els and the impact of rotation, it is useful to summarise 



the effects low metallicity has on the evolution of mas- 
sive stars. A lower metallicity implies a lower luminos- 
ity which leads to slightly smaller convective cores. This 
can be seen in Table |2 by comparing the core masses of 
non-rotating the 20 Mq models at different metallicities. 
A lower metallicity implies lower opacity and lower mass 
losses (as long as the chemical composition has not been 
changed by burning or mixing in the part of the star one 
considers). So at the start of the evolution lower metal- 
licity stars are more compact. This can be seen in the 
Hcrzsprung-Russell (HR) diagram (Fig. left) where the 
lower metallicity models have bluer tracks during the main 
sequence. They also lose less mass as can be seen by look- 
ing at the final masses in Tabled The lower metallicity 
models also have a ha rder time reaching the r ed super- 
giant (RSG) stage (see IMaeder fc Mevnetll20oH for a de- 
tailed discussion). The non-rotating model at Z — 10 -3 
becomes a RSG only during shell He-burning (see Fig. 
and the lower metallicity non-rotating models never reach 
the RSG stage. As long as the metallicity is above about 
Z = 10~ 10 , no significant differences have been found in 
non-rotating models. Below this metallicity and for metal 
free stars, the CNO cycle cannot operate at the start of 
H-burning. At the end of its formation, the star therefore 
contracts until it starts He-burning because the pp-chains 
cannot balance the effect of the gravitational force. Once 
enough carbon and oxygen are produced, the CNO cycle 
can operate and the star behaves like stars with Z > 10 -10 
for the rest of the main sequence. Shell H-burning still dif- 
fers between Z > 10~ 10 and metal free stars. Metal fre e 
stellar models are presen ted inlChieffi fe Lim ongi ll2004h . 
iHeger fc Wooslevl lho02l) and lUmeda fc Nomotol fcOO.fl . 



3.2. Rotation effects 

How does rotation change this picture? At all metallicities, 
rotation usually increases the c ore sizes, the lifetimes, the 
lumi nosity and the ma s s loss . IMaeder fc Mevnetl l|200ll) 
and IMevnet fc Maederl l)2002|) have already studied the 
impact of rota tion down to a me t allicity of Z = 10~ 5 . 
They find in IMaeder fc Mevnetl l|200l|) that rotation 
favours a redward evolution and that rotating models can 
reproduce the observed ratio of blue to red supergiants 
(B/R) in the small Magellanic cloud. We can see in Fig. 
H (middle) that the rotating models at Z = 10 -5 be- 
come RSGs during shell He-burning (See Fig. [2J. This 
does not change the ratio B/R but changes the structure 
of the star when the SN explodes. At Z = 10~ 8 (Fig. [l] 
right), the 20 Mq models do not become RSG. However 
other mass models do reach the RSG stage and the 85 Mq 
model even becomes a W R star of type WO (see below). 
IMaeder fc Mevnetl l|200lf) also find that a larger fraction of 
stars reach break-up velocities during the evolution. This 
will be further discussed in Sect. 13. 4J) . 

In IMevnet fc Maederl l|2002|) . they show that low 
metallicity (Z = 10~ 5 ) models have strong internal f2- 
gradients, which favours an important mixing. This mix- 
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Table 2. Initial parameters of the models (columns 1-5): mass, metallicity, rotation velocity [kms -1 ], total angular 
momentum [10 53 ergs] and v ini /v clit . Lifetimes [Myr] (6-8): total and core hydrogen and helium burnings. Various 
masses [Mq] (9-14): final mass, masses of the helium, carbon-oxygen, silicon and iron cores and the remnant mass. 



M ini 


Z m 




Jtol "mi/Ucrit 


T life 




THc 


M flnal 


M a 


Mco 


M Si 


M Fc 


M a 

iU rcm 


20 


2e-2 


300 


0.36 


0.44 


11.0 


10.1 


0.798 


8.7626 


8.66 


6.59 


2.25 


1.27 


2.57 


20 


le-3 


000 




0.00 


inn 


o n9 
y.uz 


n 87^ 


1 Q ^£7 


D.Oo 




1. 10 


i i « 

1.10 


9 ni 


20 


le-3 


300 


0.34 


0.39 


11.5 


10.6 


0.813 


17.1900 


8.32 


6.24 


2.27 


1.33 


2.48 


20 


le-5 


000 




0.00 


9.80 


8.86 


0.829 


19.9795 


6.24 


4.28 


1.67 


1.18 


1.98 


20 


le-5 


300 


0.27 


0.34 


11.1 


10.2 


0.806 


19.9297 


7.90 


5.68 


1.99 


1.30 


2.34 


20 


le-5 


500 


0.42 


0.57 


11.6 


10.6 


0.812 


19.5749 


7.85 


5.91 


2.18 


1.35 


2.39 


20 


le-8 


000 




0.00 


8.96 


8.24 


0.598 


19.9994 


4.43 


4.05 


1.91 


1.05 


1.92 


20 


le-8 


300 


0.18 


0.28 


9.98 


9.20 


0.610 


19.9992 


6.17 


5.18 


1.96 


1.29 


2.21 


20 


le-8 


600 


0.33 


0.55 


10.6 


9.71 


0.703 


19.9521 


4.83 


4.36 


2.01 


1.29 


2.00 


09 


le-8 


500 


0.80 


0.08 


30.5 


26.8 


3.24 


8.9995 


1.90 


1.34 






1.21 


40 


le-8 


700 


1.15 


0.55 


5.77 


5.31 


0.402 


35.7954 


13.5 


12.8 


2.56 


1.49 


4.04 


60 


le-8 


800 


2.41 


0.57 


4.55 


4.19 


0.321 


48.9747 


25.6 


24.0 






7.38 


85 


le-8 


800 


4.15 


0.53 


3.86 


3.50 


0.322 


19.8677 


19.9 


18.8 


3.19 


1.84 


5.79 


a estimated from the CO core 
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Fig. 1. HR diagrams: (left) non-rotating 20 Mq models, (middle) 20 M© models with Z = 10 5 and (right) Z = 10 
models with average rotation. 



ing leads to primary nitrogen production during He- 
burning by rotational diffusion of carbon and oxygen into 
the H-burning shell. Their results already point out that 
the primary nitrogen yields strongly depend on the ini- 
tial rotation velocity. This dependence is further studied 
in Sect. H 

iMevnet et alJ (|2006h present the evolution of 60 Mq 
models at Z = 10 -8 and confirm the effects that were pre- 
dicted in their previous papers. The fast rotating model, 
with an initial rotational velocity of 800 km s -1 reaches 
break-up, becomes a red supergiant and produces signif- 
icant amount of primary nitrogen. The model becomes a 
Wolf-Rayet (WR) star due to large mass losses during 
the RSG stage. These effects are further discussed below 
and models at Z = 10~ 8 with different initial masses are 
presented. 

3.3. Strong impact of mixing at very low metallicities 

At solar metallicity and metallicities higher than about 
Z = 10~ 5 , rotational mixing increases the helium and CO 



core masses (see Table |2J). In particular, the oxygen yield 
is increased. The i mpact of mixing on m odels at Z = 10 -8 
(and at Z = see lEkstrom et al.ll2006h is however differ- 
ent for fast rotation (ui n i =600-800 km s _1 ). The impact 
of mixing on the structure and convective zones is rep- 
resented in the Kippenhahn diagram for the Z = 10 -8 
models (see Fig. |2J. During hydrogen burning and the 
start of helium burning, the impact of mixing is the same 
as at higher metallicity: mixing increases the core sizes 
and mixing of helium above the core suppresses the inter- 
mediate convective zones linked to shell H-burning. The 
difference from higher metallicity models occurs during 
He-burning. As shown in Fig. El (left) f° r t ne 20 Mq with 
Vim =600 km s -1 , primary carbon and oxygen are mixed 
outside of the convective core into the H-burning shell. 
Once the enrichment is strong enough, the H-burning shell 
is boosted (the CNO cycle depends strongly on the carbon 
and oxygen mixing at such low initial metallicities). The 
shell then becomes convective, as can be seen in Fig. 
which is a zoom of the Kippenhahn diagram. The boost 
phase, which could look like an instantaneous event in Fig. 
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log ID (Time until core collapse) [yr] 
H- 60 , v= 800 km s 1 , t= le-08 
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i H 6 He 4 C g Ne _g Si _ 4 

log l0 (Time until core collapse) [yr] 
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Ne _ 2 Si 
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Fig. 2. Stellar structure (Kippcnhahn) diagrams, which show the evolution of the structure as a function of the time 
left until the core collapse. The initial parameters of the models are given on top of each plot. The coloured zone 
correspond to the convective zones and the symbols of the burning stages are given below the time axis. 
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M= 20 M G , v ro = 600 km s"> , Z= le-08 



s 




5.61 5.605 5.6 5.595 



log 10 (Time until core collapse) [yr] 

Fig. 3. Stellar structure (Kippenhahn) diagram of the 20 
Mq with =600kms _1 at Z = 10~ 8 . The central con- 
vective zone is created by core He-burning and the in- 
termediate convective zone is created by shell H-burning. 
This is a zoom around the time when the shell H-burning 
is boosted due to strong mixing of carbon and oxygen from 
the core. The graph shows the apparition of the convec- 
tive zone induced by shell H-burning and that the mass 
of the helium burning core is smaller after than before the 
H-burning shell boost. 



the profiles after. These profiles show that the carbon 
and oxygen brought to the shell H-burning are trans- 
formed into primary nitrogen. The bulk of primary ni- 
trogen is however produced later in the evolution. It hap- 
pens when the shell H-burning deepens in mass during 
core helium burning for the very massive models and dur- 
ing shell He-burning for the 20 Mq model. The deepen- 
ing of the convective H-burning shell is c aused by dy- 
nami cal shear instabilities (see Sect. 2.3 in iHirschi et alJ 
|2004 and refer ences therein). Dynamical shear instabil- 
ities take place just below the bottom of the convective 
zone due to the strong differential rotation between the 
convective H-burning shell and the layers below (deeper 
in the star). Dynamical shear instabilities occur on a dy- 
namical time scale as opposed to secular shear and they 
therefore induce a fast mixing just below the convective 
zone. Note that the dynamical shear instabilities are not 
influenced by mean molecular weight gradients. The con- 
sequences are the same as during the first boost. Carbon 
and oxygen are mixed inside the H-burning shell and a 
new boost occurs, producing this time more nitrogen be- 
cause the newly mixed material is richer in carbon and 
oxygen. Figure El (right) shows the abundance profiles 
after the second boost. The abundance of nitrogen does 
not change later during the advanced stages and the pre- 
SN profile for nitrogen (see Fig. 113(1 stays the same. The 
total production of primary nitrogen is further discussed 
in Sect. Rotational mixing also influences strongly the 
mass loss of very massive stars as is discussed below. 



^ (third line, right), is fully revealed in this zoom. The cal- 
culations have been repeated with different time steps to 
verify that the results did not depend on the numerical de- 
tails. The evolution of the final model was followed with 
2000 time steps between the time axis measures of 5.607 
and 5.596. 

In response to the shell boost, the core expands and 
the convective core mass decreases. At the end of He- 
burning, the CO core is less massive than in the non- 
rotating model (Fig.0 third line, left and right). The yield 
of 16 being closely correlated with the mass of the CO 
core, it is therefore reduced due to the strong mixing. At 
the same time the carbon yield is slightly increased (See 
Table It is interesting to note that the shell H-burning 
boost occurs in all the different initial mass models at 
Z = 10~ 8 as can be seen in Fig. [21 (noticeable by the 
appearance of a strong intermediate convective shell due 
to H-burning and the corresponding sharp decrease in the 
He-burning convective core mass). This means that the 
relatively "low" oxygen yields and "high" carbon yields 
are produced over a large mass range at Z = 10~ 8 . This 
could be an explanation for the possible high [C/O] ratio 
observed in the most metal poor halo stars (ratio between 
the surface abunda nces of carbon an d oxygen relative to 
solar; see Fig. 14 in lSnite et alJl200fi|) . 

Figure El (left) shows the abundance profiles before 
the shell H-burning boost and Fig. El (middle) shows 



3.4. Mass loss 

Mass loss becomes gradually unimportant as the metal- 
licity decreases in the 20 Mq models. At solar metallic- 
ity, the rotating 20 Mq model loses more than half of its 
mass, at Z = 0.001, the models lose less than 15% of their 
mass, at Z = 10 -5 less than 3% and at Z = 10 -8 less 
than 0.3% (see Table 12^1. livlevnet et all lj200rt) show that 
the situation can be very different for a 60 Mq star at 
Z = 10~ 8 . Indeed, their 60 Mq model loses about half 
of its initial mass. About ten percents of the initial mass 
is lost when the surface of the star reaches break-up ve- 
locities during the main sequence. The largest mass loss 
occurs during the red supergiant (RSG) stage due to the 
mixing of primary carbon and oxygen from the core to 
the surface through convective and rotational mixing. The 
large mass loss is due to the fact that the star crosses the 
Humphreys-Davidson limit. 

What happens in the models calculated in this study? 
First let us study the mass loss at break-up. Figure 0] 
presents the evolution of the surface velocity and of the 
ratio of the surface angular velocity to the critical angular 
velocity, fi/fi cr it f° r the 85 Mq (red solid line), 40 Mq 
(dotted black line) and 20 Mq (dashed blue line) mod- 
els with fast rotation velocities, v- ln i =600-800 km s _1 , at 
Z = 10 -8 . It shows that the 20 Mq model only reaches 
break-up velocities at the end of the main sequence (MS) 
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Fig. 4. Evolution of rotation for the 85 Mq (red solid line), 
40 M Q (dotted black line) and 20 M Q (dashed blue line) 
models with fast rotation velocities, Vim =600-800 km s _1 , 
at Z = 10~ 8 : (left) surface equatorial velocity, and (right) 
ratio of the surface angular velocity to the critical angular 
velocity, f2/f2 cr i t . 

and therefore does not lose mass due to this phenomenon. 
However, more massive models reach critical velocities 
early during the MS (the earlier the more massive the 
model). The evolution of rotation for the 60 Mq model 
is very similar to the 40 Mq model and is therefore not 
shown here for the clarity of the plot. The mass lost due 
to break-up increases with the initial mass and amounts 
to 1.1, 3.5 and 5.5 M Q for the 40, 60 and 85 Mq models 
respectively (see the top solid line in Fig. going down 
during the MS). At the end of core H-burning, the core 
contracts and the envelope expands, thus decreasing the 
surface velocity and f2/fi cr it- The mass loss rates becomes 
very low again until the star crosses the HR diagram and 
reaches the RSG stage. At this point the convective en- 
velope dredges up CNO elements to the surface increas- 
ing its overall metallicity. As said in Sect. 12.21 the total 
metallicity, Z, is used (including CNO elements) for the 
metallicity dependence of the mass loss. Therefore depend- 
ing on how much CNO is brought up to the surface, the 
mass loss can become very large again. The CNO brought 
to the surface comes from primary C and O produced in 
He-burning. As described in the above subsection, rota- 
tional and convective mixing brings these elements into 
the H-burning shell. A large fraction of the C and O is 
then transformed into primary nitrogen via the CNO cy- 
cle. Additional convective and rotational mixing is neces- 
sary to bring the primary CNO to the surface of the star. 
The whole process is complex and depends on mixing (see 
Fig. |SJ). Of particular importance is the surface convec- 
tive zone, which appears when the star becomes a RSG. 
This convective zone dredges-up the CNO to the surface. 
For a very large mass loss to occur, it is necessary that 
the star becomes a RSG in order to develop a convective 
envelope. It is also important that the extent of the con- 
vective envelope is large enough to reach the CNO rich 
layers. Finally, the star must reach the RSG stage early 
enough (before the end of core He-burning) so that there 
will be time remaining to lose mass. Figure shows the 
evolution of the effective temperature as a function of the 



M= 85 M Q , u rot = 800 km s" 1 , Z= le-08 




5.36 5.34 5.32 5.8 e 5.28 



log 10 (Time until core collapse) [yr] 

Fig. 5. Stellar structure (Kippenhahn) diagram of the 85 
Mq model at Z = 10~ 8 . This is a zoom covering the 
period when C and O are transported from the core to 
the surface and induce a strong mass loss. The central 
convective zone is related to core He-burning. The main 
intermediate convective zone is caused by shell H-burning 
and the surface convective zone appears when the star 
becomes a RSG. 

central helium mass fraction. This figure shows that the 9 
and 40 Mq models reach the RSG stage only after the end 
of helium burning, so too late for a large mass loss. The 60 
Mq model reaches the RSG stage during He-burning. It 
would therefore have time to lose large amounts of mass. 
However, the dredge-up is not strong enough. The 85 Mq 
model becomes a RSG during He-burning earlier than the 
60 Mq model. The dredge-up is stronger for this model 
and the surface CNO abundance becomes very high (see 
Fig. 03 bottom) . The series of models presented here con- 
strain the minimum initial mass for significant mass loss 
(more than half of the initial mass) to be between 60 and 

85 Mq. 

The dependence on mixing of the lower initial mass 
for a large mass loss to occur can be estimated by com- 
paring th e 60 Mp, mode l calcu lated here and the one pre- 
sented bv IMevnet et al.l ((2006) . The model calculated by 
iMevnet et al.l l)2006(l . which does not include overshooting 
and uses a different prescription fo r the horizontal diffu- 
sion coefficient, D n l|Maederl 2003), loses a large fraction 
of its mass (and becomes a WR star with high effective 
temperature) just be fore the end of cor e helium burn- 
in g (see Fig. 4 f rom IMevnet et al.l 12006^ . The D h used 
in IMevnet et al.l <|2006^> . compared to the Dh used in the 
present calculations, tends to allow a larger enrichment 
of the surface in CNO processed elements. This different 
physical ingredient explains the differences between the 
two 60 Mq models. The fact that, out of two 60 Mq mod- 
els, one model does not lose much mass and the other 
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Fig. 6. Evolution of the effective temperature as a func- 
tion of the central helium mass fraction for the Z = 10 -8 
models. Models evolve from the top right corner to the top 
left corner during the MS. They then evolve to the right 
during helium burning and downwards when they become 
RSG. The 20 M Q track resembles the 9 Af© model track 
except that it doesn't reach the RSG stage at the end of 
He-burning. 



model with a different physics just does, could mean that 
the minimum initial mass for the star to lose a large frac- 
tion of its mass is around 60 Mq. 

3.5. WR stars at very low metallicities 

The lower limit for WR star formation is roughly the same 
as for large mass loss to occur, i.e. probably about 60 Af Q , 
since models calculated here that lose a lot of mass be- 
come WRs. Therefore our models of single rotating mas- 
sive stars produce WR stars at very low metallicities. The 
85 Mq model even becomes a WO type WR star. SNe of 
type lb and Ic are therefore predicted to ensue from the 
death of single massive stars at very low metallicities. The 
next question is whether or not these stars can produce 
long and s oft gamma ray bursts (GRBs) via the collap- 
sar model Figure shows the angular 

momentum distribution at the pre-SN stage. This figure 
shows that the co re has enough angular momentum to pro- 
ducc a GRB (sec lYoon fc Langerll2005t fWooslev fc Heeerl 
120061: Iffirschi et al.ll2005al forlnore details on GRB pro- 
genitors). The effects of magnetic fields are however not 
included in the calculations. It is therefore possible that 
the core will not retain enough angular momentum to form 
an accretion disk. It is interesting to point out here that 
the outer part of the star also contains sufficient angu- 




0.4 0.6 0.S 

m/M lol , 1^=19.86773 



Fig. 7. Profile of the specific angular momentum at the 
pre-SN stage for the 85 M Q model at Z = 10~ 8 (solid 
line). The red dotted line shows the minimum angular 
momentum necessary in order to form an accretion disk 
around a rotating black hole. The blue short dashed and 
green long dashed lines show the minimum angular mo- 
mentum necessary for a maximally rotating and a non- 
rotating black hole respectively. 

lar momentum to form an accretion disk, which is less 
likely to be removed due to the effects of magnetic fields. 
In a very optimistic outcome, one could imagine that the 
model could produce two jet episodes, one when the core 
collapses and one when the outer parts collapse, since for 
such a high initial mass the central black hole could swal- 
low the whole pre-SN structure. 

4. Stellar yields of light elements and comparison 
with observations 

4.1. Stellar yields 

The stellar yields are calc ulated using the same formulae 
as in lHirschi et al.l ll2005br) . The wind contribution from a 
star of initial mass, to, to the stellar yield of an clement i 
is: 



TOB W ' nd 
'"fim 



r (m) 



M(m,t)[X?(m,t) - X?]dt 



(5) 



where r(m) is the final age of the star, M(m, t) the mass 
loss rate when the age of the star is equal to t, Xf(m,t) 
the surface abundance in mass fraction of element i and 
X® its initial mass fraction (see Table fj) . The pre-SN 
contribution from a star of initial mass, to, to the stellar 
yield of an element i is: 



pre-SN 
m Pun 



i(t) 



[Xi(m r ) - X°]dm r 



m(rcm) 



(6) 
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Fig. 8. Abundance profiles for the 40 (top), 60 (middle) and 85 (bottom) M Q models. The pre-SN and wind (yellow 
shaded area) chemical compositions are separated by a red dashed line located at the pre-SN total mass (-Mfi nal ), 
given below each plot. 



Raphael Hirschi: Pre-SN evolution and primary nitrogen @ very low Z. 



11 



Table 3. Initial mass, metallicity and velocity and total stellar yields" (mp^f - N +irip"™ ld ) of the models. All masses 
and yields are in solar mass units and velocities are in kms -1 . For the metallicity, Z ini is in loglO units (-3 corresponds 
to Zi n i = 10 -3 ) except for Z = 0.02 (©). Z is the total metal content and is def ined by: Z= 1 — Ai H — A"3 Ho — A<i He . 
These are the yields to be used for chemical evolution models using Eq. 2 from llMaederlllQQ^ . 



A f 

iVl 




Wini 


rie 


10 4 


12 ri 


Li 


14 m 
IN 


iSn 


i7n 


U 


INe 


z 


20 





300 


1.62 


-3.36 


0.433 


1.01e-3 


4.33e-2 


2.57 


-2.75e-6 


-1.96e-4 


4.26e-2 


3.98 


20 
20 


CO CO 


000 
300 


2.47 
2.11 


-2.84 
-3.31 


0.373 
0.676 


2.58e-5 
2.84e-5 


3.31e-3 
3.10e-3 


1.46 
2.70 


-5.48e-7 
4.83e-7 


-1.12e-5 
-1.89e-5 


2.48e-3 
1.04e-2 


2.38 
3.77 


20 
20 
20 


-5 
-5 
-5 


000 
300 
500 


2.50 
2.34 
2.26 


-3.32 
-3.91 
-4.08 


0.370 
0.481 
0.648 


1.93e-7 
2.42e-6 
1.53e-5 


4.27e-5 
1.51e-4 
5.31e-4 


1.50 
2.37 
2.59 


3.05e-7 
3.40e-7 
4.79e-7 


-9.43e-8 
5.27e-7 
5.49e-6 


2.11e-5 
2.74e-3 
1.07e-2 


2.30 
3.35 
3.54 


20 
20 
20 


00 QO 00 


000 
300 
600 


2.27 
2.03 
3.15 


-4.17 
-4.37 
-4.50 


0.262 
0.381 
0.823 


2.27e-4 
1.80e-6 
5.55e-3 


8.52e-3 
1.20e-4 
5.90e-2 


1.20 
1.96 
1.35 


1.94e-7 
1.70e-8 
1.73e-5 


-2.15e-0 
2.14e-7 
2.52e-7 


1.85e-6 
5.48e-5 
7.72e-5 


2.14 
2.97 
2.49 


09 
40 
60 

85 


00 00 00 00 


500 
700 
800 
800 


1.43 
6.01 
8.97 
16.8 


-1.76 
-9.10 
-13.3 
-2.00 


0.082 
1.79 
3.58 
7.89 


1.35e-4 
6.31e-2 
5.00e-4 
5.60e-l 


2.53e-3 
1.87e-l 
4.14e-2 
1.75e+0 


5.85 
5.94 
12.8 
12.3 


5.64e-7 
6.31e-5 
6.26e-6 
6.66e-4 


4.07e-5 
4.74e-7 
3.56e-7 
4.95e-5 


1.86e-4 
1.51e-4 
1.66e-3 
1.55e-3 


0.143 
9.64 
17.1 
25.5 



a Note that the corresponding ejected masses can be calculated by adding the initial composition given in Table multiplied 
by the mass interval, the mass boundaries of which (initial and remnant masses) are given in Table |5] 



Table 4. Initial mass and velocity and stellar wind yields (mpf^ d ) of the Z = 10 8 models. All masses and yields 
are in solar mass units and velocities are in kms -1 . 



M 


"ini 


4 He 


3 He 


12 C 


13 C 


i4 N 


ib 


i7 


18q 


^Ne 


Z 


09 


500 


2.80e-5 


-8.34e-9 


5.86e-08 


1.72e-09 


2.53e-08 


2.33e-08 


6.19e-12 


1.20e-ll 


9.77e-12 


1.10e-7 


20 


600 


2.36e-4 


-1.04e-6 


-3.27e-ll 


5.06e-13 


3.08e-10 


-2.54e-10 


5.30e-13 


-7.22e-13 


2.72e-16 


7.59e-9 


40 


700 


3.29e-l 


-1.03e-4 


5.34e-03 


8.06e-04 


3.63e-03 


2.42e-03 


1.05e-06 


2.18e-09 


2.35e-07 


1.17e-2 


60 


800 


1.21e+0 


-2.72e-4 


1.80e-05 


4.94e-06 


6.87e-04 


5.48e-05 


7.94e-08 


-9.49e-ll 


4.80e-08 


7.25e-4 


85 


800 


2.00e+l 


-1.64e-3 


6.34e+00 


5.60e-01 


1.75e+00 


3.02e+00 


6.66e-04 


4.95e-05 


1.46e-03 


1.16e+l 



a Note that the corresponding ejected masses can be calculated by adding the initial composition given in Table Q multiplied by 
the mass interval, the mass boundaries of which (initial and final masses) are given in Table |2] For heavy elements with yields 
larger than le-6, ejected masses and stellar yields are essentially the same here since the initial total metallicity is le-8. 



where m(rem) is the remnant mass estimated from the CO 
core mass (see lMaederlll992|) . m(r) the final stellar mass, 
Xf the initial abundance in mass fraction of element i 
and Xi(m r ) the final abundance in mass fraction at the 
lagrangian mass coordinate, m r . The pre-SN contribution 
is calculated at the end of Si-burning. Therefore the con- 
tribution from explosive nucleosynthesis is not included. 
However, elements lighter than n eon are marginally modi- 
fied by explosive nucleo synthesis ijChiefii fc Limongl2003l 
iThielemann et alJl99r^ and are mainly determined by the 
hydrostatic evolution. The total stellar yield of an element 
i from a star of initial mass, m, is then : 



mp. 



pre-SN 



(7) 



m Plm corresponds to the amount of element i newly syn- 
thesised and ejected by a star of initial mass m. One notes 
that, using the above formulae, negative yields are ob- 
tained if the star destroys an element (see yields of 3 He). 
The total stellar yields for the chemical elements which are 
not significantly affected by the evolution beyond the cal- 
culations done in this work are presented in Table |3| The 
stellar wind contribution is presented for the Z = 10 -8 
models in Tabled The SN contribution to the yields (lay- 
ers of the stars ejected during the SN) is not presented but 



the pre-SN abundance profiles are presented in Figs.[S]and 
1131 The 60 Mq model was evolved until neon burning and 
the 9 Mq model until carbon burning. This means that, 
for these two models only, the abundance profiles can still 
vary in the central regions before the collapse. Since the 
remnant mass probably contains a large fraction of the 
CO core for these two models, the SN contribution to the 
yields should not be affected by these variations. 

The 20 Mq models are expected to produce neutron 
stars and eject most of their envelope. On the other hand, 
stars more massive than about 40 Mq on the ZAMS are 
expected to form black holes directly a nd not to eject any - 
thing during or after their collapse ijHeger et al.l l2003|) . 
The present 40 and 60 Mq models probably follow this 
scenario. If this is the case, the wind contribution is the 
only contribution to be taken into account. The outcome 
is uncertain for the 85 Mq model because the final mass 
is only about 20 Mq but the alpha and CO core masses 
are very large (see Table 0). This model could produce 
a GRB, in which case jets would be produced and some 
iron rich matter would be ejected. In Fig. |5| the two pos- 
sible outcomes for very massive stars are compared. On 
the left, the stellar yields include the SN contribution and 
on the right, the stellar yields only include the wind con- 
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Initial mass [M ] Initial mass [M Q ] 



Fig. 9. Total stellar yields as a function of the initial mass 
of the models with Z = 10 -8 and (left) assuming the SN 
ejects matter for M > 40 Mq and (right) assuming direct 
black hole formation and no matter ejection by the SN for 
M > 40 M Q . 

tribution for stars with M > 40 Mq. As can be expected, 
the two outcomes give very different yields for the 40 and 
60 Mq since they do not lose much mass before the SN 
explosion. For stars above about 60 Mq, the difference is 
much smaller because the strong winds peal off most of 
the CNO rich layers before the final collapse. 

4.1.1. Metallicity dependence of the CNO yields 

The metallicity dependence of the yields can be studied 
using the 20 M P) models. In deed when the yields are con- 
volved by a ISalpeterl ill 9 Sift like initial mass function, 20 
Mq is not far from the average massive star. Note that 
it is however uncertain whether a standard IMF applies 
at very low metallicities. Recent studies tend to argue 
that the IMF is top heavy but very massive stars (VMS, 
M > 140M Q ) are not necessary to reproduce observations 
(see for example ISchneider et alJ l2006t iTumlinson et alJ 
12004 - 

The yields for 12 C, 14 N and 16 are presented in Fig. 
1101 The most stringent observational constraint at very 
low Z is a very high primary 14 N production. This requires 
extremely high primary 14 N production in massive stars, 
of the order of 0.1 M fl per star (~0.15 used in the 
heuristic model of IChiappini et ai1l2005|) . In Fig. we 
can see that only the model at Z — 10 -8 and with a fast 
rotation (u ini =600 km s -1 ) gets close to such high values. 
The bulk of 14 N is produced in the convective zone created 
by shell hydrogen burning (see Sect. YS.'d} . If this convec- 
tive zone deepens enough to engulf carbon (and oxygen) 
rich layers, then significant amounts of primary 14 N can 
be produced (0.01^0.1 Mq). This occurs in both the non- 
rotating model and the fast rotating model but for differ- 
ent reasons. In the non-rotating model, it occurs due to 
structure rearrangements similar to t he third dredge-up 
at th e end of carbon burning (see also IChieffi fc Limongil 
2004, for a similar process in Z=0 models). In the model 
with v- ln i =600 km s -1 it occurs during shell helium burn- 
ing because of the strong mixing of carbon and oxygen 
into the hydrogen shell burning zone. Another interesting 



feature is the possible upturn of the [C/O] ratio observed 
at very low metallicities (ratio between the surface abun- 
dances of carbon and oxygen relative to solar; see Fig. 14 
ISpite et al.1l2005(l . Indeed, the models at Z = 10~ 8 with a 
fast rotation have high C and low O yields compared to 
the Z = 10 -5 models and could reproduce such an upturn. 
As explained in Sect. IH-.HI the reason is the reduction of 
the He-burning core mass after the shell H-burning boost 
induced by strong mixing. As can be seen in Fig. [21 this 
occurs for all initial masses calculated in the Z = 10~ 8 
series. The stellar yields of the fast rotating models calcu- 
lated here were used in a galactic chemical evolution model 
and suc cessfully reproduce the early evolution of CNO el- 
ements l|Chiappini et alJ l20061. This is a good argument 
in favo ur of fast rotation at v ery low metallicities. Note 
that in IChiappini et all (|2006), the SN contribution was 
included in the stellar yields for all masses. The impact 
of including only the wind contribution for masses above 
40 Mq will be studied in the future. The impact of the 
important primary nitrogen production and of the other 
yields on the initial composition and therefore on the evo- 
lution of the next stellar generations and their yields is an 
interesting aspect that will also be studied in the future. 

4.1.2. Z = 1(T 8 models 

The stellar yields of the fast rotating Z = 10~ 8 mod- 
els are presented in Fig. It shows that the signifi- 
cant (above 0.01 Mq) production of primary nitrogen 
occurs for the entire mass range. For massive stars with 
40 < M < 60 Mq, the yields depend on whether or not 
the SN contributes to the total yields as discussed earlier. 
The production of nitrogen is accompanied by a produc- 
tion of 13 C and 17 (and to a lesser extent 18 0). The ratio 
12 C/ 13 C is very different between the wind and the SN 
contributions. It is around 5 for the wind and more than 
100 for the SN and SN+wind contributions. This, with 
the N/C and N/O ratios are good tests to differentiate 
between the two contributions. Part of the primary nitro- 
gen also captures two ct-particles and becomes 22 Ne. The 
primary 22 Ne yields are of the order of 10 _4 M Q . 22 Ne is 
one of the main neutron sources for the weak s-process in 
massive stars. With a low initial iron content due to the 
low initial overall metallicity, the s-process could occur 
with a high neutron to seed ratio and produce surprising 
results. This will be the subject of a future study. 

4.2. Carbon rich EMP stars 

The zoo of the extremely metal p oor stars has 
been classified by iBeers fc Christliebl l|2005l) . Carbon 
rich extremely metal poor stars (CEMPs also 
called CRUMPs at a recent meeting at Tegernsee, 
http : / / www. mpa-garching . mpg . de/ ~crumps0 5) are dif- 
feren t from normal EMP star and are rarer l|Rvan et alJ 
2005). About three quarter of the CEMPs show a stan- 
dard s-process enrichment pointing to the fact that they 
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Fig. 10. 12 C (left), 14 N (centre) and 16 (right) total stellar yields as a function of the initial metallicity of the 20 
Mq models. The red solid curves corresponds to the models with fast rotation, the blue dashed curve to the very 
low metallicity models with Ui n i=300kms _1 and the dotted black cu rve to the non-rotating models. In the graph 
of the 14 N yields, the yields used in the chemical evolution models of IChiapnini et a"fl |2005) in order to reproduce 
observations is marked with the symbol CMB05. For the fast rotating models, one sees that the carbon yield increases 
slightly between Z = 10~ 5 and Z = 1CP 8 . At the same time, the oxygen yield decreases towards lower Z by a factor 
about two. 



have ac creted matter from an AGB companion in a binary 
system l)Suda et alJl2004|) . The other CEMP stars show a 
weak s-process enhancement. Their peculiar abundances 
are therefore thought to originate from the previous gen- 
eration of stars. The two most metal poor stars known to 
date, HE1327-232 6 llFrebel et all 120051 lAoki et al.ll2006t 
[ Frebel et all 120061) and HE 0107-5240 ijchristlieb et all 
l2004t iBessell et all Eool are both CEMP stars w ith a 
weak s-proces s enric hment. IChieffi fc Lim T^il ll2002|) and 
iNomoto et, all l|2005h have studied the enrichment due to 
PopIII SNe. By using one or a few SNe and using a very 
large mass cut, they are able to reproduc e the abundance 
of m ost elements 1 Limongi et all 120031 llwamoto et all 
l2005h . However they are not able to reproduce the 
nitrogen surface abundance of HE 1327-2326 without 
rotational mixing or the oxygen surface abundance of HE 
0107-5240 without mixing and fall-back mimicking an 
aspherical explosion. In this work, the impact of rotation 
is explored. HE1327-2326 is characterised by very high N, 
C and O abundances, high Na, Mg and Al abundances, a 
weak s-process enrichment and depleted lithium. The star 
is not evolved so has not had time to bring self-produced 
CNO elements to its surface and is most likely a subgiant 
(Korn, presentation at the CRUMPS meeting). A lot 
of the features of this star are similar to the properties 
of the stellar winds of very metal poor rotating stars 
l|Mevnet et alJl2006j) . HE1327-2326 could therefore have 
formed from gas, which was mainly enriched by stellar 
winds of rotating very low metallicity stars. In this 
scenario, a first generation of stars (PopIII) pollutes the 
interstellar medium to very low metallicities (TFe/Hl~-6). 
Then a PopII.5 star l|Brommll2005l lHirschill2005l iKarlssonl 
l2006h like the 40 Mq model calculated here pollutes 
(mainly through its wind) the interstellar medium out 
of which HE 1327-2326 forms. This would mean that 



HE1327-2326 is a third generation star. In this scenario, 
the CNO abundances are well reproduced, in particular 
that of nitrog en, which according the new values for a 
subgiant from iFrebel et alJ l)2006l) is 0.9 dcx higher in 
[X/Fe] than oxygen. This is shown in Fig. 1111 where the 
new abundances are represented by the red stars and the 
best fit is obtained by diluting the composition of the 
wind of the 40 Mq model by a factor 600. On the right 
side of Fig. ^2 one sees that when the SN contribution is 
added, the [X/Fe] ratio is usually lower for nitrogen than 
for oxygen. The lithium depletion cannot be explained 
by rotational and convective mixing in the massive star 
if the wind material is diluted in the ISM by a factor 
600 (600 parts of ISM for 1 part of wind material) as is 
suggested above. However, if the wind material is less 
enriched in CNO elements, a lower dilution factor would 
be necessary to reproduce the observations. Also if the 
massive star is born with a higher iron content, a lower 
dilution factor is necessary. If this dilution factor is of 
the order of unity, it becomes possible to explain the 
lithium depletion by internal mixing in the massive star. 
To investigate this possibility, more models have to be 
calculated with different initial metallicities. Although 
the existence of a lower limit for the minimum metallicity 
Z for low mass stars to form is still under debate, It is 
interesting to note that the very high CNO yields of the 
40 Mq stars brings the total metallicity Z above the limit 
for low mass star formation obtained in iBromm fc LoefJ 
l|2003h . 

For HE 0107-5240 (red hexagons in Fig.[nj, rotation 
does not help since neither the wind contribution nor the 
total contribution produce such large overproduction of 
carbon compared to nitrogen and oxygen. Possible ori- 
gins for this star are p resented in llwamoto et al.1 l|2005h 
and i Suda et all l2004|) . For the other carbon rich stars 
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Fig. 11. Composition in [X/Fe] of the stellar wind (left) 
and the mixture of the wind and SN ejecta (right) 
for the Z = 10~ 8 models. The lines represent predic- 
tions from the models. The following symbols correspond 
to the values obtained from measurements at the sur- 
face of giant CEM P stars: red hexagon s for HE 0107- 
5240, [Fe/H]~ -5.3 l|Christlieb et alJl2Q04h: areen thin tri- 
angles CS 22949-03 7 j [Fe/Hl~ -4.0 jNorris et aD l2002t 

"* " 200411 : r.msses CS 
|200^ 



iDepagne et all 120021 : llsraelian et alJ 
29498-043, [Fe/H]~ -3.5 llAoki etal 



Two values 

are given for oxygen for the last two stars. They come 
from the determinations us i ng O triplet or [OI] f orbid- 
den lines llAoki et alJl2004 llsraelian et al.lf2004ll . The 
empty t riangles JPlez fc Cohenll2005h . [Fe/H]~ -4.0 and 
squares (|Frebel et al.ll2005|) . [Fe/H]- -5.4 correspond to 
non-evolved CEMP stars. The ne w (3D/NLTE correc ted) 
estimates for HE1327-2326 from Iprebel et ahl l|2006j) are 
represented by the red stars. On the right, wind+SN ejecta 
are compared to the normal EMP st ars (green hatched 
area with lines go ing down to the right ICavrel et al.l l2004: 
ISnite et alJl2005h for the 20 Mq model and again to the 
CEMP stars for the more massive models. For this pur- 
pose, the value [O/Fe] is chosen to fall in the middle of 
the observed range (20 M©: [O/Fe]=0.7 and M > 20 M Q : 
[0/Fe]=2.8). For HE1327-2326, the best fit for the CNO 
elements is obtained by diluting the composition of the 
wind of the 40 Mq model by a factor 600. 



presented in Fig. ^] the oxygen abun dances are either 
not d etermined or still quite uncertain l)lzotov fc Thuanl 
120041). The C and N surface abundances of G77-61 
l(Plez fc Cohenll2005h could originate from material sim- 
ilar to the wind of the 8 5 Mp, . The C and N surface abun- 
dances of CS 22949-037 llNorris et al.l2002HDepagne et all 
l2002t llsraelian et al.ll2004b^ resemble the wind composi- 
tion of the 60 Mq model although some oxygen needs 
be ejected from the supernova. The enrichments in C. N 
and O are very s imilar for CS 29498-043 (around +2, see 
lAoki et alJl2004|) and a partial ejection due to the super- 
nova is necessary to explain the oxygen enrichment since in 
the winds the oxygen is usually under produced compared 
to C and N. It will be interesting to follow the evolution 
of Na, Mg and Al since the high yields of 22 Ne seem to 
indicate that there could an overproduction of these ele- 
ments in the wind (see Table and the Z = 10~ 5 model 



presented in lMevnet et al.l2006h . Since 22 Ne is also a neu- 
tron source, s-process calculations are also planned. 

5. Conclusion 

Two series of models were computed. The first series 
consists of 20 Mq models with varying initial metal- 
licity (solar down to Z = 10~ 8 ) and rotation (uj n i = 
— 600 km s -1 ). The second one consists of models with 
an initial metallicity of Z = 10~ 8 , masses between 9 
and 85 Mq and fast initial rotation velocities. The re- 
sults presented confirm the crucial role of rotation in stel- 
lar evolution and its impact in very low metallicity stars 
[|Mevnet et alJl2006(l . The evolution of the models with 
Z = 10~ 8 ([Fe/H]~ —6.6) is very interesting. In the course 
of helium burning, carbon and oxygen are mixed into the 
hydrogen burning shell. This boosts the importance of the 
shell and causes a reduction of the CO core mass. Later 
in the evolution, the hydrogen shell deepens and produces 
large amount of primary nitrogen. For the most massive 
models (M > 60 Mq), significant mass loss occurs during 
the red supergiant stage assuming that CNO elements are 
important contributors to mass loss. This mass loss is due 
to the surface enrichment in CNO elements via rotational 
and convective mixing. The models predict the produc- 
tion of WR stars for an initial mass higher than 60 Mq at 
Z = 10~ 8 and the 85 Mq model becomes a WO. Therefore 
SNe of type lb and Ic are predicted from single massive 
stars at these low metallicities. The 85 Mq model retains 
enough angular momentum to produce a GRB but the 
calculations did not include the effects of magnetic fields. 

The stellar yields are presented for light elements. 
These yields were used in a galactic chemical evolution 
model and suc cessfully reproduce the early evolution of 
CNO elements ( Chiappini et al.ll2006h . A scenario is also 
proposed to explain the abundances of the most meta l 
poor star known to date HE1327-2326 l|Frebel et al.l2005|) . 
In this scenario, a first generation of stars (PopIII) pollutes 
the interstellar medium to very low meta llicities ([Fe/H]— - 
6). Then a PopII.5 star l|Bromml liooll) like the 40 Mq 
model calculated in this study pollutes only with its wind 
the interstellar medium out of which HE1327-2326 forms. 

There are still many questions or issues that could not 
be treated in this work. It is necessary to determine over 
which metallicity range the large primary production and 
the other specific features of the Z = 10 -8 models occur. 
The impact of the yields on the initial composition and 
therefore on the evolution of the next stellar generations 
also needs to be studied. It will also be interesting to fol- 
low the evolution of Na, Mg and Al since the high yields of 
22 Ne seem to indicate that there could an overproduction 
of these elements in the wind (see Tab le 01 and the Z = 
10~ 5 model presented in lMevnet et al.ll2006|) . Since 22 Ne 
is also a neutron source, s-process calc ulations are also 
planned. The effects of magnetic fields llYoon fc Langerl 
l2005UWooslev Hegerll2006UMa,eder Rr. Mevnetll2005h on 
the results will be studied in the near future. The depen- 
dence of the mass loss rates on the metallicity, especially 
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in the RSG stage need to be further studied to see how 
the results of lvan Lo'onl (|2005l mass loss in the RSG phase 
independent of metallicity) can be extrapolated to very 
low metallicities. 
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Fig. 12. Abundance profiles for the 20 M with =600 km s^ 1 at Z — 10~ 8 : (left) before shell H-burning boost 
(induced by mixing of carbon and oxygen from the core), (middle) just after the boost and (right) after shell H-burning 
has deepened during shell He-burning. 
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Raphael Hirschi: Pre-SN evolution and primary nitrogen @ very low Z. 
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Fig. 13. Final abundance profiles. The initial parameters of the models are given on top of each plot. 



